Asymmetry dependence of proton correlations. 
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A dispersive optical model analysis of p+ 40 Ca and p+ 48 Ca interactions has been carried out. The 
real and imaginary potentials have been constrained from fits to elastic scattering data, reaction 
cross sections, and level properties of valence hole states deduced from (e, e'p) data. The surface 
imaginary potential was found to be larger overall and the gap in this potential on either side of 
the Fermi energy was found to be smaller for the neutron-rich p+ 48 Ca system. These results imply 
that protons with energies near the Fermi surface experience larger correlations with increasing 
asymmetry. 
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In the independent-particle model, nucleons in the nu- 
cleus move in a mean-field potential generated by the 
other nucleons. All nucleon levels up until the Fermi 
energy (Ep) are fully occupied, while those above are 
empty. Although this model enables an understanding 
of various aspects of nuclear structure, a full descrip- 
tion of nuclei and nuclear matter requires consideration 
of the correlations between the nucleons. These include 
short-range, central and tensor interactions and longer 
range correlations associated with low-lying collective ex- 
citations |2(. As a result, for closed-shell nuclei, single- 
particle (sp) levels below E F have an occupancy of only 
70-80% and the levels at higher energy have a nonzero 
occupancy Q. The strength of the sp levels are spread 
over energy, with narrow peaks or broad distributions 
(depending on their separation from Ep). In addition, 
there is strength at very high momentum |^. 

Although there are numerous studies of the effect of 
correlations on the properties of sp levels for nuclei near 
stability, there are only a few studies for very neutron 
or proton-rich nuclei. From neutron knock-out reactions, 
Gade et al. [5( infer the occupancy of the 0c?5/ 2 neutron 
hole state in the proton-rich 32 Ar nucleus is considerably 
reduced relative to those for stable nuclei. 

An alternate method to study sp strength is through 
the use of the dispersive optical model (DOM) developed 
by Mahaux and Sartor [f|. This description employs the 
Kramers-Kronig dispersion relation that links the imag- 
inary and real parts of the nucleon self-energy Q • This 
procedure links optical-model (OM) analyses of reaction 
data at positive energies to structural information at neg- 
ative energies. In the present work, the properties of 
proton levels in Ca nuclei as a function of asymmetry 
5= n ~a Z are investigated with the DOM. Previously mea- 
sured elastic-scattering and reaction-cross-section data 
for protons on 40 Ca and 48 Ca as well as level proper- 
ties of hole states in these nuclei, inferred from (e, e'p) 
reactions, were simultaneously fit. The dependence on 5 
is extracted and used to predict level properties of 60 Ca. 

In the DOM, the complex energy-dependent potential 
felt by the protons is comprised of a real V, volume W v 



and surface W s imaginary components, plus spin-orbit 
V so and Coulomb V c potentials, 

U{r, E) = -V(r, E) + V so (r) + V c (r) 

-iW v (E) f(r,r v ,a v )+Aia s W s {E)^-f(r,r s ,a s ). (1) 

Wood-Saxon form factors f(r, r%, ai) = [1 + e a » ] _1 
are used. The real part of the nuclear potential is as- 
sumed to be given by two terms 

AV{r,E) = V HF (E) f(r,r HF ,a HF ) + AV(r,E) (2) 

where Vhf has a smooth energy dependence arising from 
the nonlocality or momentum-dependence of the micro- 
scopic self-energy. The dispersive correction AV has vol- 
ume and surface parts, 



AV(r, E) = AV v (E)f(r, r v ,a v )+4a s AV s (E)—f(r, r s ,a s ), 

ar 

(3) 

and is related to the imaginary potential through the 
dispersion relationship, i.e., 



AVi(E) = ^ J Wi{E') 



1 



E F 



dE' 



(4) 

where i — v, s and P stands for the principal value. The 
dispersive corrections are a result of coupling to non- 
elastic channels. The surface term accounts for the influ- 
ence of low-lying collective states and giant resonances. 

The form for the imaginary potential must take into 
account the dominance of surface and volume absorp- 
tion at low and high positive energies, respectively. In 
addition around E F the imaginary potentials must be 
zero. The potential should be approximately symmetric 
around Ep, but further away it must become asymmetric 
as there are a finite number of hole states. In this work 
we assume 

Wv {E) = Av {E-Erf + B* + AWnA ' L {E) (5) 
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FIG. 1: (Color online) Calculated and experimental elastic- 
scattering angular distributions of the differential cross sec- 
tion da/dfl and analyzing power. Left panels shows results for 
p+ 40 Ca, right for p+ 48 Ca. Each curve is labeled by the pro- 
ton energy in MeV. For display purposes, successively higher 
energy curves and data for da/dQ are scaled down by an addi- 
tional factor of 4.5. For the lowest energies, compound-elastic 
contributions were also included in the fits. 



where the energy- asymmetric correction AWnm.{E) i s 
derived from nuclear-matter considerations [6( . The sur- 
face potential was taken as the difference of two functions 
that cancel at large energies, i.e., 

AW S (E) = oj(E, A l a ,Bl,c, 0) - uj(E, A%B 2 s ,c, Q), (6) 

X 4 



u(E,A,B,c,Q)=AG(X) 



,-cX 



X 4 + B 4 



(7) 



where 9 (X) is Heavyside's step function, X — \E — Ep\ — 
Q, A 2 S = A\e- C Q, and Q = B l a + AB. 

The Hartree-Fock potential is often assumed to de- 
crease linearly or exponentially with energy. We took 
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FIG. 2: (Color online) Total reaction cross sections are dis- 
played as a function of proton energy. 



the form 



Vhf(E) 



2A 



HF 



l + exp(£W (E-E f )/A H f) 



(8) 



which is approximately linear around Ep and becomes 
more exponential at larger energies. This form provided 
a reasonable location for the 0si/ 2 level in 40 Ca. 

The parameters of the DOM were fit for both 40 Ca 
and 48 Ca from a large set of data covering both positive 
and negative proton energies. For 40 Ca, 14 experimental 
elastic-scattering angular distributions for energies from 
18 to 135 MeV [3, HE El E and seven data sets 
for the analyzin g p ower measured at energies from 21 
to 80 MeV [H Hi [H H [H EH were included. For 
48 Ca the fitted data included 14 angular distributions 
and seven sets of analyzing power at energies from 8 to 
65 MeV [II [II [II |ij . Reaction cross sections for both 
targets were taken from the tabulations of Bauhoff [2l| . 
For the 0d$/2, l s i/2j an< i 0g?3/2 proton holes states, the 
mean level energies, r.m.s. radii, spectroscopic factors, 
and widths inferred from measured (e, e'p) cross sections 
[22I E3| were also included. Lastly, the fits considered 
the mean ener gies of the Of 7/2 and O/5/2 (for 48 Ca only) 
particle levels |24| . 

The seven geometric parameters 

(thf, a>HF=av,r v , r s ,a s ,r so ,a so ) defining the form 
factors were varied, but kept identical for both targets. 
Similarly the fit parameters B 2 , AB and c defining 
the decay of the W S (E) at large energies, as well the 
magnitude of the spin-orbit term, were also set identical 
for both targets. Only Ahf, Bhf, A],Bl, A v , and B v 
were allowed to differ for each isotope. 

The final fit to the experimental data is shown in 
Figs. lll4l and the fitted potentials are displayed in Fig.|SJi. 
It was found possible to obtain similar quality fits with 
different potentials, however a robust feature of all good 
fits was that the magnitude of the surface imaginary po- 
tential ( A],) was larger and its width B 1 (of the minimum 
around Ep) was narrower for the neutron- rich Ca nu- 
cleus. There was however an ambiguity in determining 
the rate at which this potential diminished at large en- 
ergies. This ambiguity is coupled with an ambiguity in 
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FIG. 3: (Color online) Fitted level properties for the 0d 5/ / 2 , 
lsx/2, and 0d :i / 2 proton hole states. The left (a,c,e) and right 
(b,d,f) panels display the fits for 40 Ca and 48 Ca, respectively. 
The fitted quantities include a,b) the root-mean-squared ra- 
dius Rr.m.s, c,d) the spectroscopic factor S expressed as a 
percent of the independent-particle-model value, and e,f) the 
widths F of these levels. 



determining the magnitude and the rate of increase of 
W V (E). 

Tornow et al. 2±g fitted similar data for 40 Ca with the 
DOM using a W s (E) which decreased slowly and was 
still substantial at the highest energy considered in the 
work. With such a slow diminishing of W s (E) , one can 
also obtain good fits to the 48 Ca data, however the fitted 
W v (E) potentials are substantially different for 40 Ca and 
48 Ca. On the other hand, if W s (E) is made to diminish 
faster, these differences can be reduced to zero. 

Ambiguities in determining potentials in standard OM 
fits are well known, however volume integrals of the po- 
tentials have been shown to be better defined Com- 
parisons of the imaginary volume integral, Jw(E) = 
J W(r, E)dr, obtained from the OM fits in the referenced 
experimental studies, indicate that J\y is larger below 
E ~50 MeV in 48 Ca. However for higher energies, there 
is no discernible difference for the two isotopes. Thus 
for 48 Ca, if its larger surface potential is still significant 
in this higher energy region, W v {E) must be smaller to 
produce similar values of Jw- We believe this result is ar- 
tificial, and therefore a solution where W s (E) diminishes 
faster is preferable. In any case, it is not possible to con- 
strain any difference in W v (E) for the two nuclei. Thus, 
for the final fit we have taken W v (E) to be independent 
of asymmetry. This is consistent with global OM fits 
|2(i l27j which have a significant asymmetry dependence 
for W s (E), but none for W v (E). Theoretically some 
asymmetry dependence of W v (E) would be expected and 



FIG. 4: (Color online) Comparison of experimental proton 
single-particle levels iHllllIllilliilMlIllto D ° M ° al " 
culations. For the levels indicated with the solid dots, their 
energies were included in the fits. The dashed lines indicate 
the Fermi energy. 



higher-energy 8 Ca data would provide sensitivity to this. 

Calculated and experimental sp level schemes for 40 Ca 
and 48 Ca are displayed in Fig. 0] Apart from the levels 
included in the fit (indicated with the solid dots), the 
other known levels are well reproduced. The 0si/ 2 level 
of 40 Ca is very wide and even though the DOM prediction 
for its energy is l ow, it lies within the experimentally 
determined width |28j . 

Present DOM calculations have been extrapolated to 
60 Ca assuming the parameters Ajjf, Bhf, A\, and 
B] vary linearly with S. The extrapolated energy- 
dependence of W S (E) is shown in Fig. as the thin- 
lined curve and the predicted sp level scheme is displayed 
in Fig. The surface dispersive correction is large for 
this nucleus and its effects are quite apparent. The lev- 
els in the immediate vicinity of Ep are focused closer to 
Ep, increasing the density of sp levels. A reduced gap 
between the particle and hole valence levels implies that 
the closed-shell nature of this nucleus has diminished and 
proton pairing may be important. The levels further from 
Ep have been pushed away and as a result there are big 
gaps between the 0pi/ 2 and the 0d 5 / 2 and also between 
the Fpi/2 and Ogg/ 2 levels. 

The occupation probabilities as defined in Ref. [|| are 
plotted in Fig. [5J;. The results for 48 Ca agree with the 
theoretical work of Ref. ■ For proton hole states just 
below Ep, the occupation probabilities have decreased 
for the more neutron-rich 48 Ca while the opposite is true 
for the particle states. This trend is further accentuated 
in our extrapolation to 60 Ca. On the other hand, the 
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FIG. 5: (Color online) Quantities derived from DOM fits, a) 
The energy dependence of the Hartree-Fock potential Vhf, 
and the surface W s and volume W v imaginary potentials, b) 
The radial dependence of the effective mass at Ef- c) The 
occupation probabilities (points indicate sp levels). For all 
plots, the thick (back), medium (red) and thin (green) curves 
give the results for 40 Ca, 48 Ca, and 60 Ca, respectively. 

more deeply-bound 0s 1 / 2 , 0p 3 / 2 , and Opi/2 levels show 
very little sensitivity. Their occupancies are more sensi- 
tive to the volume imaginary component whose asymme- 
try dependence was not constrained. 

Our observation that the occupancies of valence pro- 
ton hole states are reduced in neutron-rich 48 Ca can be 
compared to the reduced occupancy of the valence neu- 
tron hole state in the proton-rich 32 Ar inferred by Gade 
et al. 5]. Thus a preponderance of one type of parti- 
cle, reduces the occupancies of valence hole states for the 
other type. This indicates that correlations are stronger 
for these valence nucleons. Recent calculations of asym- 
metric nuclear matter by Frick et al. [3j| also predict 
such a result as a purely volume effect. Its origin is 
probably similar in both cases reflecting that p-n inter- 
actions are stronger than n-n or p-p, partly because of 
the tensor force. Thus, protons in neutron-rich systems 
are more strongly correlated as illustrated from the in- 
ferred spectroscopic factors of 65%, 56% and 50% for the 
0^3/2 proton level in 40 Ca, 48 Ca and 60 Ca, respectively. 
Conversely for neutrons, the opposite is expected. 

The nucleon effective mass to* (r, E)/m = 1 — 
dV(r, E)/dE (to is the nucleon mass) at E=Ep inferred 
from this work is displayed in Fig. 03. Only the surface 
contribution around r—A-5 fm has been constrained and 
it increases significantly with asymmetry. This suggests 
an asymmetry dependence of the surface component of 
the level-density parameter. 

In conclusion, the properties of proton single-particle 



states in the vicinity of Ep for °Ca and Ca have been 
studied with a comprehensive dispersive-optical-model 
analysis of elastic-scattering and bound-level data. The 
analysis indicates that the imaginary surface potential is 
~50% larger and the minimum around the Fermi energy 
is narrower for the neutron-rich 48 Ca nucleus. This im- 
plies that, with increasing asymmetry, the occupancies of 
proton levels vary more smoothly across the Fermi sur- 
face, a consequence of increased correlations. The present 
observations and those of Gade et al. || can be under- 
stood from the larger strength of p-n relative to the p-p 
and n-n interactions. Hence, protons (neutrons) experi- 
ence larger (weaker) correlations in neutron-rich matter. 
The reversed is true for proton-rich matter. 

This work was supported by the U.S. Department of 
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